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Abstract. A system of Schwinger-Dyson equations for pseudoscalar four-dimensional 
Yukawa model in the two-particle approximation is investigated. The simplest it- 
erative solution of the system corresponds to the mean-field approximation (or, 
equivalently, to the leading order of 1/A-expansion) and includes a non-physical 
Landau pole in deep-Euclidean region for the pseudoscalar propagator A. It is ar- 
gued, however, that a full solution may be free from non-physical singularities and 
has the self-consistent asymptotic behavior p 2 A ~ C log~ 4 ^ 5 jfe. An approximate 
solution confirms the positivity of C and the absence of Landau pole. 



PACS number: 11.10. Jj. 

1 Introduction 



A definition of asymptotic behavior at large momenta for the strictly renormalized 
four-dimensional models of quantum field theory (QFT) up to now is an unsolved 
problem. The unique exception is the theory with non-Abelian gauge interaction, 
which is asymptotically free in the framework of renormalization-group-improved 
perturbation theory. A solution of the problem of asymptotic behavior for other 
models requires going out the framework of the coupling-constant perturbative ex- 
pansion. The first attempt to define the asymptotic behavior in QFT was made 
by Landau and coworkers in the 1950s. This investigation was based on the ap- 
proximate solution of Dyson equations with the summation of leading logarithms, 
and the result for quantum electrodynamics was as follows: the photon propagator 
included the non-physical singularity in Euclidean region of momenta [1]. Then 
similar singularities were indicated in the model with Yukawa interaction [2] . Such 
singularities in the Euclidean region violate general principles of QFT and are a 
serious problem for these models. Further development has demonstrated that 
these non-physical singularities arise practically inevitably in the framework of any 
known non-perturbative methods: at the renormalization-group summation, in the 
frameworks of 1/iV-expansion and mean-field expansion, etc. 2 

A widespread opinion is formulated as a triviality of the quantum field mod- 
els that is not asymptotically free in the sense of the improved coupling-constant 
perturbative expansion. There is a rigorous theorem [4] that the four-dimensional 
scalar field theory with 4 interaction on the lattice does not have an interacting 
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continuum theory as its limit for zero lattice spacing, i.e. the theory is trivial. How- 
ever, this argument is not fully conclusive due to an uncertainty of the continuous 
limit in this model [5]. 

In spite of the serious evidence for the triviality of lattice scalar theory, the 
situation with triviality today is as vague as before, and recent papers in this topic 
maintain incompatible statements. While mainstream works confirm the triviality 
scenario, Suslov in a series of works(see [6] and references therein) argues the non- 
trivial behavior for </> 4 -theory and quantum electrodynamics in the strong-coupling 
region. In any case, a study of the triviality problem requires a non-perturbative 
tool. 

In this paper we investigate a new non-perturbative approximation for a model 
with Yukawa interaction - two-particle approximation, or 2PA. This approximation 
was proposed for scalar fields in [7] , and it is the first non-trivial step of a sequence 
of general n-particle approximations, which tends to the exact infinite system of 
Schwinger-Dyson equations (SDEs) at n — > oo. 

The structure of the paper is as follows: in section 2, the necessary notations 
and definitions are given; SDEs for the generating functional of Green functions are 
introduced in the formalism of a bifocal fermion source 3 . We consider using of the 
bifocal source as a convenient choice of the functional variable. In particular, this 
variable is very convenient for the construction of the mean-field (MF) expansion, 
which is presented in section 2. The existence of the Landau pole in the boson 
propagator of the leading approximation of the MF expansion is also demonstrated 
in this section. 

In section 3, a general construction of the approximation scheme for the system 
of SDEs is given. The renormalization of the system of equations is made, iterative 
solutions and some supplement simplifications are discussed. 

In section 4, the asymptotic solution of the system at large Euclidean momenta 
is presented and the asymptotic behavior of the propagators at large momenta is 
discussed. The boson propagator in this model possesses self-consistent behavior. 
Conclusions are presented in section 5. 

2 Schwinger-Dyson equations, mean-field approx- 
imation and Landau pole 

Consider the theory of a Dirac fermion field ip interacting with a pseudoscalar boson 
field in a four-dimensional (1+3) space with the Lagrangian 

C = tp(id - m)V - ^</>(M 2 + d 2 )4> + gi>Tijxj> (1) 

3 A formalism of the bilocal source was first elaborated in QFT by Dahmen and Jona- 
Lasinio [8]. 
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Here T = i-y 5 and d = j^d^. The renormalizability in all orders of the coupling- 
constant expansion requires to supplement Lagrangian (1) with an additional term 
A</> 4 which corresponds to the self-interaction of the boson field. This term en- 
sures a renormalization of the boson-boson scattering amplitude. In this work we 
are not concerned with this amplitude, and the renormalization of approximations 
considered below does not require including the corresponding counter-term. For 
this reason we do not include the quartic interaction into consideration. Therefore, 
we shall consider the restricted Yukawa model neglecting the self-interaction of the 
boson field, and the obtained results should be treated as the first step to the study 
of the asymptotic behavior in a realistic model of the boson-fermion interaction. In 
other words, we shall consider the case A = only. Also we shall believe m > 
and M 2 > and, therefore, do not discuss in this work the problem of dynamical 
mass generation by Yukawa interaction (see [9]). 

The generating functional of Green functions can be written as a functional 
integral 



G = D(ip, ip, 4>) expi 



dx C 



dxdyip(y)ri(y,x)ip(x) + / dxj(x)<j>(x) 



, (2) 



where j{x) is a single boson source and rj{x,y) is a bilocal fcrmion source. 

The translational invariance of the functional integration measure leads to the 
functional-differential SDEs for generating functional G. In terms of the logarithm 
Z = i log G these equations are: 



S(x — y) + (id x - m)i 



5Z 



Sr)(y,x) 



S 2 Z 



+ i- 



5Z 



5Z 



5r)(y,x)5j(x) 5-q{y,x) 5j(x) 

SZ 



I 



dxir)(x, X\) 



SZ 
5j(x) 



dxi A c (x - xi)j(xi) - gA c (x - xi)tr( T — 

V 6r}(xi, 



8Z 



Xi) 



Here A c = (M 2 + d ) . We define also the fermion propagator 



S(x -y) = i 



SZ 



the boson propagator 



A(x -y) = TT 



Sr](y,x) 
S 2 Z 



5j(y)Sj(x) 

the two-particle (four-point) fermion function 

S 2 Z 



v=j=o 



7 i x y 

22 [ x' y> 



5r](y',x')5r](y,x) 



(3) 
(4) 

(5) 
(6) 

(7) 



n=3=o 
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and the three-point function 



G 3 {z\x,y) = i- 



S 2 Z 



(8) 



r)=j=0 



5t](y,x)Sj(z 

Differentiations of SDE (4) over ij and j give us the SDE for the three-point function 

G 3 (z\x, y) = -gj d Zl A c (z - Zl )T Z 2 ( * » ) (9) 

and the SDE for the boson propagator 

A(x-y) = A c {x-y)+ig J d yi ti[T G 3 (x\y u Vl )] A c ( Vl - y). (10) 

Excluding with the help of SDE (4) a differentiation over j in SDE (3), we obtain 
at j = the SDE for the generating functional: 



+^A c (x-a;i)r 



S(x - y) + (id x - m)i 
SZ 



5Z 



-tr r 



Sr)(y,x) 
SZ 



+ 



dxi^ irj(x,xi) 
5 



SZ 



-tr r- 



5r](y,xi) 
SZ 



+ 



5-q{y,x) \ 5r](x 1 ,x 1 )J 5r)(y,x) V 6r}(xi,xi) 



)]} <"> 



which contains only the derivatives over the bilocal source r\. Switching off the 
source r\ in (11), we have the SDE for the fermion propagator 

(m - id x ) S(x -y) = S(x - y) + ig 2 j d Xl A c (x - Xl )TZ 2 ^ ^ ^ ^ Y (12) 

A differentiation of (11) over r\ gives us (with the source being switched off) the 
SDE for the two-particle fermion function 

(m - id x ) Z 2 f * V \+5(x- ?/) S(x' - 



V) = 



ig- 



' Jdx 1 {(rs(x- y ))A c (x-x 1 )z 2 (^ x x ) *) ^t + ta c (x- Xi )z 3 ^ xi xi jr}(13) 



Here Z 3 = i^f | _ is the three-particle (six-point) fermion function. The deriva- 
tion of equation (13) implies that tr(j 5 S) = 0, i.e. we suppose parity conservation. 

To construct the MF expansion, we consider as a leading approximation for 
equation (11) the equation 



S(x - y) + (id x - m)i 



SZ 



MF 



i J dxi^rj(x,Xi) 



SZ 



MF 



+ g 2 A c {x-x 1 )T 



SZ 



MF 



-tr 



8r)(y,x) 

SZ MF 



Sr](y,x) L Sr](xi,xi) 



} (14) 
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We call (14) the MF approximation, since this equation gives us the same equations 
for propagators and the two-particle funcion as the MF expansion for the generating 
functional (see [7] and references therein). 
The MF fermion propagator is 

S = S C , (15) 

where S c = (m — id)^ 1 . 

Equation (14) gives the equation for the two-particle function 



3>( * V y ,^+S c {x~y')S c {x'-y) = 
= ig 2 J d Xl dx 2 (S c (x - xi) T S c { Xl - y)) A c ( Xl - x 2 )Z 2 ( J X *^j T 



(16) 



whose solution is 



x y 
x' y' 



-- -S c (x - y') S c (x' - y) + 
+ j d Xl dx 2 (S c (x - Xl ) T S c ( Xl - y))f MF (xi - x 2 )(S c (x' - x 2 ) V S c (x 2 - y')), (17) 



where in the momentum space 

S«7 4;V|tlW (18) 

and 

LV) = / ^^ c ( P + q )TS c { q )T] (19) 

is the single fermion loop. 

Taking into account SDEs (9) and (10), we obtain the MF boson propagator 

A(p 2 ) = ^ f MF (p 2 ). (20) 

The above formulae contain divergent integrals and should be renormalized. 

The simplest method of renormalization of the above equations is the direct 
application of a regularization procedure in the spirit of Bogolyubov R-operation 
[10]. The unrenormalized mass operator is 

*(p 2 ) = A- V) - A" V) - -W 2 L C (P 2 ). (21) 

Then the renormalized mass operator is defined as 

a r (p 2 )^ r ega(p 2 ) = - l g 2 L^p 2 ), (22) 
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where 

L%(j?) = L V) - L c (0) p 2 (L c )>(0) = ^£dz log[l - z(l z)^}. (23) 

is the rcnormalized fermion loop ( for the easement of the following calculations we 
choose the normalization point at zero momenta). As a result, the renormalized 
boson propagator A r is 

A; 1 (p 2 ) = M 2 -p 2 -ig 2 L c r (p 2 ) (24) 

We can also renormalize the MF approximation by introducing counter-terms 
in the Lagrangian. In correspondence with the standard recipe, we consider (1) 
as the renormalized Lagrangian, where ip, <j>, m, M and g are now the renormalized 
fields, masses and coupling, and add counter-terms 

A£ = -^fVV-^V (25) 

which absorb the divergences. 

The full Lagrangian C b = C + AC can be written as 

C b = C + AC = Mid - m b )ip b - -4> h {Ml + 2 )(f> b + gbMMb (26) 

where 

tp b =ip, (f> b = y/z^(t), g b = -y=, m b = to, M 2 = . (27) 

Then all the above calculations are reproduced with Lagrangian (26), and the nor- 
malization conditions are imposed on the renormalized propagator A r . 

The normalization conditions for the propagator A r (p 2 ) = z^ 1 A b (p 2 ) are 

A- 1 (0) = M 2 , A A - 1 | p2=0 = l. (28) 

These conditions define the mass-renormalization counter-term SM 2 and the field- 
renormalization constant z$. Then the renormalized boson propagator is defined 
by equation (24) as above. 

As it follows from equation (24), the renormalized boson propagator A r pos- 
sesses a non-physical singularity (Landau pole) in the Euclidean region p 2 < at 
the point p 2 = —p 2 = M 2 , where M 2 is a solution of the equation 

M 2 + M 2 -ig 2 L c r (-M 2 ) =0. 

This equation has a solution at any positive g 2 . As was yet noted in the introduc- 
tion, the same Landau pole arises in the calculations of the renormalized amplitude 
by other methods: in the frameworks of 1/A-expansion and renormalization-group 
summation. 
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3 The system of SDEs and two-particle approxi- 
mation 



The system of fermion SDEs is an infinite set of equations for n-particle fcrmion 
functions Z n = iS n Z/Srj n \, n= Q. The first SDE is equation (12). The second SDE 
is equation (13). The nth SDE is the (n - l)th derivative of SDE (11) with the 
source being switched off and includes a set of functions from one-particle fermion 
function S to (n + l)-particle fermion function Z n+ \. 

We call " the n-particle approximation of the system of SDEs" the system of n 
SDEs in which the first n — 1 equations are exact and the nth SDE is truncated 
by omitting the (n + l)-particle function. It is evident that the sequence of such 
approximations goes to the exact set of SDEs at n — > oo. The one-particle ap- 
proximation is simply equation (12) without Z 2 . This approximation has a trivial 
solution which is a free propagator. The two-particle approximation is a system of 
equation (12) and equation (13) without Zy. 



(m - id x ) Z 2 ( X x , ^ ^+6(x- y') S(x' - y) = 
ig 2 J dx 1 {TS{x-y))A c {x-x 1 )Z 2 (^ X J ) T 



(29) 



which includes S and two-particle function Z 2 . This nonlinear system will be the 
object of the present investigation. 

The idea of these approximation scheme is very simple and natural. However, 
the calculations became more and more complicated at each following stage: e.g., 
the three-particle approximation is a system of three nonlinear equations for the 
propagator, the two-particle function and the three-particle function. 

Another view to the origin of system (12) and (29) is based on a modification 
of the MF expansion of section 2 with taking into account a particular solution 
of functional-derivative equation (11). It is easy to see that the SDE (11) has the 
simple solution 

ZM = y i Z 2 -v 2 -*s- v , (30) 

where functions Z 2 ( , , ) and S(x — y) satisfy the system of equations (12) 

and (29). To be an exact solution of the functional-derivative equation (11), this 
system should be supplemented by one more nonlinear equation for Z 2 . The system 
of three equations for the two functions S and Z 2 are overfull and seemingly has 
not physically meaningful solutions. However, the third equation does not play a 
role for the construction of the modified MF expansion. The construction of this 
expansion for scalar field theory can be found in work [7]. This construction can 
be generalized also to fermion fields. At the nth step of this expansion, we have a 
closed system of linear integral equations, and therefore this scheme is much less 
complicated in the calculational sense in comparison to the above scheme of the 
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n-particle approximations. Equations (12) and (29) are the basic approximation 
for this expansion. 

Equations (12) and (29) are the system of nonlinear equations for the functions 
S and Z 2 . In equation (29), the two-particle function Z 2 can be considered as a 
functional of S, and the "solution" of this equation can be easily found: 

l I ) = S c {x-y')S{x'-y) 

+ J dx x dx 2 (S c (x - Xl )T S( Xl - y))f(xi - x 2 ){S{x' - x 2 )T S c (x 2 - ?/)) (31) 

Here 

f(x -y) = -ig 2 A c (x - y) + ig 2 J dx 1 dx 2 A c (x - x 1 )L(x 1 - x 2 )f(x 2 - y) (32) 

and L(x) = tr[S c (x) T S(~x) T] is the fermion loop operator. 
In momentum space: 

1 =^A-\ P *) + L(p% (33) 



f(p 2 ) 9 

L(p2) = Swf tr[sc{p +q)T s{q) r] - (34) 

Taking into account equations (31)-(34), (9) and (10), we obtain for the boson 
propagator A the following equation in momentum space: 

A-V) = M 2 -p 2 -ig 2 L{p 2 ). (35) 

From equation (12) and equations (31)-(35) we have the equation for the fermion 
propagator 

S- 1 (p) = m-p + ig 2 K(p), (36) 

where 

K(p) = J ^TS c (p-q)TA( q ). (37) 

The system of equations (34)-(37) is the system of unrenormalizcd SDEs in the 
2PA. 

The renormalization of equations (35) and (36) can be performed in correspon- 
dence with the general recipe of section 2. 

If we define the unrenormalized mass operators as 



Z(p) = S- 1 (p)-S- 1 ( P ) =ig 2 K(p) , 
a(p 2 ) = A-V) - A-V) = -W 2 L( P 2 ) 
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then the renormalized mass operators are 



E r (p) = rcg E(p) = ».g 2 A-(p) - (0) - p 



dK 



p=0 



^ <j r (p 2 ) = reg ^(p 2 ) = -i 9 2 (l{ P 2 ) - i(0) -p 2 L'(0)) 
and the system of renormalized equations for propagators are 

Sr 1 (p) = m - p + ig 2 K r (p) 
A-V) = M 2 - p 2 - ig 2 L r {p 2 ) 



(39) 



(40) 



where 



K r (p)=K(p)-K(0)-p^ 



p=0 



(41) 
(42) 

The renormalization with counter-terms can be performed by introducing the 
counter-term Lagrangian 



L r (p 2 ) = L(p 2 )-L(0)-p 2 L'(0). 



A£ = — l)tp idip — 8m ipip — 



z<p 



l60 2 6- SM ^ 2 



2 2 

and full Lagrangian (26). The bare quantities in (26) are now 



^b = y/z^tp, <Pb 



9b 
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8m 



Zip \J Z<fi 



Mi 



Z-0 



M 2 + SM 2 



(43) 



(44) 



The bare mass operators Ef, and a b in the 2PA are given by formulae (38) with 
substitutions E — > E(,, a — > <Jb, g — > gb, m — > m^, M 2 
conditions (28) for the propagator A r (p 2 ) and for the fcrmion propagator 



M 2 , etc. Normalization 



S- 1 (p = Q) = 



m. 



dS-\p) 



dp 



= -1 



(45) 



p=0 



define the counter-terms, and the system of renormalized equations will be system 
(40) again. 

Note that an iteration of equation for A r in (40) with = S c leads to the 
MF propagator (24). So the MF approximation and the equivalent leading-order 
l/A^-expansion are contained in the 2PA as the first iteration. 

Equations (40) are the system of nonlinear integral equations for the propaga- 
tors. The most interesting problem is to look for the asymptotic behavior of the 
solution of system (40) at large Euclidean momenta. In the large-momenta region, 
an essential technical simplification is possible, namely, one can replace in integrals 
(34) and (37) the function S c by a massless function —l/p: 



Scip-q) 



_ f A F(q) 
J (2tt) 4 p 



(46) 
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Then it is possible to use the well-known formula 



f(q 2 e) 



(2tt)4 (p - q )l 16tt 2 



/<: i-oo 

f(q 2 e)qldq 2 e + f{ql)dql 
hi 



(47) 



(48) 



This massless-integration approximation (46) is quite usual in investigations in the 
deep-Euclidean region, though rigorous arguments for its validity can be done for 
the asymptotically-free models only [5]. In the general case, this approximation 
should be considered as a plausible conjecture, which needs further investigations. 
Formula (47) highly enables the calculations and, as a major point, permits us to 
go from integral equations to differential ones (see below). 

Equations (40) in the massless-integration approximation in the Euclidean re- 
gion are 

' <V\) = 1 - I"" d <le(Pl - & A r(q 2 e ) 

A-\pl) = M 2 +p 2 e -^ jf dqt{pl q 2 ) 2 Jgfr 
Here a is defined by S^ 1 with formula 

S- 1 =6- op, 

and b = m in the massless-integration approximation. Euclidean normalization 
conditions are 

A^O) = M 2 , (A-!)'(0) = 1, a(0) - 1. (49) 
Introducing the dimensionless quantities 

o2 



M 2 ' M 2 ' w P 2 e A r 

system (48) can be written as follows 



(50) 



M 2 +tia 2 (ti) 



Let us discuss the iterative solutions of system (50). Due to the nonlinearity 
of equations an essential moment is a way of iterations. Consider two schemes of 
iterations. The hrst scheme is a choice a(°) = 1 as a leading order, and then the 
calculation of from the second equation with a — a^°\ As was pointed above, 
this scheme corresponds to the MF approximation and leads to the Landau pole in 
the boson propagator A r at point t^. At this point h^it^) = 0, and the following 
calculation of from the hrst equation is possible only at t < t^. For this reason 
the subsequent calculations in the framework of this scheme became problematical 
in the region of large t. 
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Another iterative scheme is a choice of the leading order as ho = 1 + j , and the 
following calculation of a\ from the first equation with h — h a : 

< "-^( 1 + l) Sl °* +( > + 1 + l& + ^ < 51 > 

Likewise h^ in the first scheme, function ai(t) equals to zero at some point to, but 
consequences of this fact are quite different. Since /i 2 + ta\ > at any t > 0, 4 then 
h\ can be calculated for any t, where at t — > oo 

hi ~ 4 log log £ > 0, (52) 

i.e., hi has a self-consistent asymptotic behavior. Certainly, a finite number of 
iteration cannot define the actual asymptotic behavior, but the given considerations 
demonstrate the obvious preference of the second scheme in comparison with the 
usual MF approximation. This preference consists in the absence of the Landau pole 
in the fermion propagator, which is defined by a\ of equation (51) at ^ 2 > 0. This 
circumstance is essentially used for the construction of an approximate solution in 
the following section. 



4 Asymptotic behavior 

System (50) (multiplied to t 2 ) after threefold differentiations is reduced to the 
system of differential equations 



diS\ b a ) ifyP th ,-■_), 

d3 (t 2 h) - 3 2 » ( ' 

dt5\ L n ) — 47r 2 ^ +a i t 



At t — > system of integral equations (50) gives us: 

«M = l-^+^ 2 + ^) (54) 

and 

h { t) = \ + l-^t + { t 2 ) (55) 

These formulae give us boundary conditions for system (53) at the point t = 0. 
At large t, system (53) has the asymptotic solution 5 : 

a ~ A log 1/5 1 , h ~ B log 4/5 1, AB = - . (56) 

327T 



A case /i =0 needs a special consideration and is not discussed here. 
5 Note that system (53) has exact solution a = Ct~ 1/2 , h = , C 2 = — f ^t 2 . This 

imaginary solution is not, of course, a solution of integral equations (50). 
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Differential equations (53) do not fix the signs of coefficients A and B. These 
signs have the principal meaning and define the physical situation, described by 
2PA. If A > (B < 0), then function h(t) applies to zero at some point. This case 
corresponds to the presence of Landau- type singularity in the boson propagator, i.e., 
the situation is similar to the above physically unsatisfactory MF approximation. 
Oppositely, if A < 0(B > 0) , then the situation corresponds to self-consistent 
asymptotic behavior of the boson propagator. In this case, function a(t) has a zero 
at some point, but it does not lead to a Landau pole in the propagator for massive 
fermions. 

System (53) of nonlinear differential equations is rather difficult for detailed 
analytical investigation, and below we make some approximations and natural sim- 
plified suppositions, which will enable the construction of an approximate solution 
and the definition of the coefficients A and B. 

Firstly note, that function v — ah has quite definite asymptotic behavior at 
t — > oo: 

5<7 2 

v = ah~-^\ogt. (57) 

Since v ~ j at t — > 0, this function changes the sign in some point to and v(to) = 0. 
We shall suppose the uniqueness of this point. Consider the behavior of the solution 
near the point to- After the change of variable x — logi, the first equation of system 
(53) can be written as 

Here a' = Suppose in a vicinity of zero point x = logio : 

v w v'(x )(x - x ). (59) 

Note that supposition (59) is fulfilled for iterative solution (51), which will be used 
in the construction of the approximate solution (see below). 

Then, going to the variable z = x — xo = log j-, we obtain for a the linear 
differential equation 

a" 1 + 3a" + 2a' = k-, (60) 

Z 

where k — — , / , . Since v decreases, k > 0. With the substitution 

a = e- z u (61) 
and by changing the variable £ = j, we obtain for u the Meijer equation [11] 

Ktt!4-«!-w!-3>«!i-°- 

whose solution is 

u = dm + C 2 u 2 + C3M3, (62) 
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where 



Ul 



-^(-tI Vtlf )•— 



z 

T 



1 - k/2 
|;0 



w 3 = G 



2 

T 



1 

0,1; i 



fc/2 n 



(63) 



At z — > we have Ui ~ z 2 , U2 ~ z, M3 ~ z°. Consequently, if C3 7^ 0, function a 
docs not change sign at point x , and A > 0. If C3 = 0, C2 7^ 0, function a changes 
sign, and A < 0. This case corresponds to the above-mentioned situation of self- 
consistent asymptotic behavior. The case C3 = C 2 = corresponds to touching for 
a and a pole singularity for h. 

A definition of coefficients Cj needs some boundary conditions. These boundary 
conditions should be connected with boundary conditions at t = and can be 
defined on the basis of some approximate solution in the region of small t. The 
linearized version type of equation (60) is tightly connected with the asymptotic 
behavior at large t and apparently cannot be applied in the region of small t. For 
this reason, in the pre- asymptotic region we shall use iterative solution (51), which 
satisfies supposition (59) and boundary conditions (54). In other words, we believe 
a = ai at t < to (x < xq). In the region of large t, we use asymptotic formula (57) 

for v and believe v = —^^r{ x — x a) at t > t (x > x ). Correspondingly, we have 
for a at large t equation (60), whose solution is given by equations (61)-(63) with 
k = 2/5. 

As the boundary conditions we shall use the conditions of "a smooth join" at 
point x — x (z = 0): 

a(z = 0) = a(x ) = ai (so) = 0, a' z (0) = a[(xo), <4(0) = a'{(x ). (64) 



From these boundary conditions, we define 



r(i.2) 



{a , ;{x ) + 2a' 1 {x )), C 2 



T(0.7) 



a'^xo), C 3 = 0. (65) 



Here T(k) is the Euler gamma- function. Simple calculation shows that both C\ 
and C2 are negative. Asymptotic expansions of Meijer functions (see [11]) give us 
the behavior of solution at z — > +00: 



= Ciui + C 2 u 2 ~ e z (|) V5 ^= (C 2 + d) 



(66) 



and, consequently, the asymptotic behavior of this approximate solution is given 
by formulae (56) with 

A = — S t= (Ci + C 2 ) < (67) 



and B > 0. At t > 1 



2 6 / 5 0F 



-4^-0.77^, S«3.24, 



(68) 



13 



and 



0.308 



Pi log 



4/5 vl 



(69) 



5 Conclusions 

Our results demonstrate that two- particle approximation of the system of Schwinger- 
Dyson equations essentially differs in the asymptotic deep-Euclidean region of mo- 
menta in comparison with the MF approximation. Instead of physically unsatis- 
factory behavior with the Landau pole in the Euclidean region, which occurs for 
the MF approximation, or for the leading term of 1/iV-expansion, the boson prop- 
agator in the two-particle approximation has self-consistent asymptotic behavior, 
which is similar to the asymptotically free behavior. Certainly, inclusion of the 
scalar self-action into consideration can notably vary the results in the asymptotic 
region. From this point of view, the obtained results should be considered as the 
first step of investigation of the full model, which will include Yukawa interaction 
and the self-action of the scalar field. 
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